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ABSTRACT 

The variation of the specific intensity across the stellar disc is essential input 
parameter in surface brightness reconstruction techniques such as Doppler imaging, 
where the relative intensity contributions of different surface elements are important 
in detecting starspots. We use phoenix and ATLAS model atmospheres to model 
lightcurves derived from high precision (S/N ~ 5000) HST data of the eclipsing binary 
SV Cam (F9V + K4V), where the variation of specific intensity across the stellar disc 
will determine the contact points of the binary system lightcurve. For the first time 
we use x 2 comparison fits to the first derivative profiles to determine the best-fitting 
model atmosphere. We show the wavelength dependence of the limb darkening and 
that the first derivative profile is sensitive to the limb-darkening profile very close 
to the limb of the primary star. It is concluded that there is only a marginal differ- 
ence (< lex) between the x 2 comparison fits of the two model atmospheres to the HST 
lightcurve at all wavelengths. The usefulness of the second derivative of the light-curve 
for measuring the sharpness of the primary's limb is investigated, but we find that the 
data are too noisy to permit a quantitative analysis. 
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1 INTRODUCTION 

Limb darkening effects in stellar atmospheres have impor- 
tant implications throughout stellar astrophysics where a de- 
termination of the surface brightness distribution is impor- 
tant. Recent work using Doppler imaging and micro- lensing 
events has shown that commonly-used analytical limb dark- 
ening la ws fail to match stellar observations at the limb of 
the star llThurl et al.l2004l:lBarnes et all2004) . Other micro- 
lensing results (lFieldsll2003ft show that the intensity predic- 
tions from model atmospheres are discrepant in the case of 
a K-giant. 

Surface brightness reconstruction techniques such as 
Doppler imaging and eclipse mapping rely on the informa- 
tion content of surface areas with diff erent distan c es fro m 
the rotation axis (see review bv ICollier Cameron! (|2001)). 
To detect starspots Doppler imaging uses the relative in- 
tensity contributions, calculated from model atmospheres, 
of the different surface elements. To reconstruct an accu- 
rate surface brightness distribution it is essential to know 



how parameters, such as the limb darkening, can alter the 
intensity values across the stellar disc. 

High inclination eclipsing binary systems can be used as 
probes to determine the variation of specific intensity at the 
stellar limb. If limb darkening showed a smooth transition in 
specific intensity at the limb of the star, the contact points 
of eclipses would appear less abrupt and slightly displaced in 
phase relative to models with limb darkening laws derived 
from plane-parallel atmospheres, where the cutoff is very 
sharp. The sharp cutoff in plane parallel atmospheres results 
from the optical depth of the rays being infinite at the limb. 

In November 2001 we were awarded 9 orbits of 
HST/STIS time to eclipse-map the inner face of the F9V 
primary of the totally eclipsing binary SV Cam. SV Cam 
(F9V + K4V) is a synchronously rotating RS CVn binary 
with a period of 0.59d. We obtained spectrophotometric 
lightcurves of 3 primary eclipses with a signa l-to-noise ra- 
tio of 5000. The first analysis of these data, bv l.Teffers et alJ 
(2005), determined the radii of the primary and secondary 
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stars. When the resulting lightcurve was subtracted from the 
observed data, the re sidual lightcurve sh owed strong peaks 
at phases of contact. I Jeffers et alJ (1200.4 ) then showed that 
these mismatches are reduced significantly, but not elimi- 
nated, when a polar cap and a reduction in the photospheric 
temperature, to synthesise high spot coverage, are imposed 
on the image. 

As there is a significant temperature difference between 
the primary and secondary stars, the secondary star acts as 
a dark occulting disc as it eclipses the primary star. The 
variation of brightness as a function of phase reflects the 
degree of limb darkening on the primary star. In this paper 
we determine the best fitting model atmosphere by fitting 
the models to the brightness variations as the secondary 
scans the inner face of the primary star, using the first and 
second derivatives of the HST in 10 wavelength bands. We 
discuss the implications these results have on results from 
Doppler Imaging. 



2 MODEL ATMOSPHERES 

In this paper, two well established stellar atmosphere 
codes are used; the phoenix model atmosphere code 
(Hauschildt et al Jl999t) that uses sp herical atmosp heres and 
the ATLAS model atmosphere code <lKurucz|ll994l) that uses 
plane parallel atmospheres. 

2.1 ATLAS 

We use the plane-parallel ATLAS9 model atmospheres from 
the Kurucz CD-ROMS (lKurucdll994l) . We integrate the in- 
tensity values over the wavelength range of our observations 
2900A to 5700A. We use temperature models from 3500 K 
to 6500 K, with 250 K interval, across 17 limb angles. The 
limb angle fi is defined by fi= cos 8, where 9 is the angle 
between the line of sight and the normal vector of the local 
surface element. The treatment of convection is based on 
the mixing length th eory with approximate overshoot (ref. 
ICastelli et alJ |l997)) with a mixing length to scale height 
ratio of 1.25. The variation of specific intensity, i.e. at /i=l,as 
a function of wavelength and limb angle is shown in Figure^ 

2.2 PHOENIX 

The general input physics set -up of the phoenix mod el at- 
mosphere code is discussed in iHauschildt et alJ 1^999). The 
main advantage of using this code is that it is based on 
spherical geometry (spherical radiative transfer) LTE rather 
than traditional plane-parallel structure. NLTE effects are 
considered to be insignificant in this application. 

The synthetic spectra are based on an extension of 
the grid of PHOE NIX model atmospheres described by 
lAllard et alJ (|2000). This extended grid includes surface 
gravities larger than log(g) > 3.0 needed f or main sequence 
stars. These models are as described by IHauschildt et all 
( 1999), but include an updated molecular line list. The mod- 
els are computed in spherical geometry with full atomic and 
molecular line blanketing using solar elemental abundances. 
In these models, the stellar mass is 0.5 Mq and the con- 
vection treatment assumes a mixing-length to pressure scale 
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Figure 1. The predicted variation in specific intensity across the 
stellar disc as a function of wavelength for the longest (5596A) 
and the shortest (3198A) wavelengths in this analysis (T=6000K, 
log g=4.5). 

height ratio of 2 with no overshooting. There are 117 syn- 
thetic spectra in total. The effective temperature runs from 
2700K to 6500 K in 100K steps at three surface gravities: 
Iog(<?) = 4.0, 4.5, and 5.0. The wavelength resolution of these 
synthetic spectra is lA. 

The variation of specific intensity as a function of limb 
angle is shown in Figure Q The difference between the 
phoenix and the ATLAS model atmospheres at the limb of 
the star results primarily from the effects of spherical geom- 
etry of the phoenix model. In models using spherical geom- 
etry, there is a finite optical depth for rays close to the limb, 
while with plane parallel models the optical depth of the rays 
is infinite at the limb providing significant intensities down 
to /i=0. Other contributing effects include overshooting and 
the mixing length ratio. 

In addition to the limb effects, Fig.Qshows that ATLAS 
intensity profiles are overall brighter than phoenix inten- 
sity profiles. We attribute this difference to the use of over- 
shooting in the ATLAS models. ATLAS models with overshoot- 
ing have been shown to have weaker limb darkening in the 
blue relative to models without overshooting, for instance, 
overshooting models have been shown to better fit solar 
limb darkening obse rvations than non-overshooting models 
ICastelli et alJ il997tl . A similar resu lt holds for Procyon (F5 
IV- V): ref. lAufdenberg et alJ l|2005h for a detailed study of 
the effects of 1-D and 3-D convection on limb darkening. 

2.3 Model lightcurves 

We u se the eclipse-mapping code DoTS llCollier Cameronl 
1997) to model the primary eclipse lightcurve for each model 
atmosphere, and to determine which model atmosphere best 
fits the complete HST lightcurve. The input data to each 
model comprises: the variation of the specific intensity as a 
function of limb angle, as shown in Figure for both model 
atmospheres considered; the primary and secondary radii; a 
reduced primary photospheric temperature and polar spot. 
Gravit y darkening is also inclu ded according to the descrip- 
tion of lCollier Cameronl Jl997f) . 

We include a reduced photospheric temperature to ac- 
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Visit 


Obs. Date 


UT 


UT 


Exposure 


No of 






Start 


End 


Time(s) 


Frames 


1 


01 November 2001 


20:55:56 


01:00:17 


30 


165 


2 


03 November 2001 


14:34:29 


18:38:21 


30 


165 


3 


05 November 2001 


09:49:17 


13:52:22 


30 


165 



Table 1. HST Observations of SV Cam, where the exposure time 
is per frame. 
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Figure 2. The comprising sections of the primary eclipse from 
each of the 3 HST visits. The flux density has been summed 
between 2900 A and 5700 A. An off-set of 4xl0 6 ergs/s/cm 2 has 
been added for clarity. 



count for a stellar surface that is peppered with small spots 
which are too small to be resolved through ec lipse mapping. 
Following the results of Jef fers e t al. (2005) we include a 
polar cap in the modelled lightcurve to optimise the fit of 
the lightcurve to the data. The method for modelling the 
photometric lightcurve to include a reduced photospheric 
temperature and a polar cap is described in Appendix A for 
the ATLAS model atmosphere. The binary system parame- 
ters are summarised in Table |!|] The lightcurve solutions for 
the two model atmospheres are in good agreement. 
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Figure 3. The variation of the HST lightcurve with wavelength. 
Each lightcurve has been plotted with an offset of 5 X 10~ 12 from 
the lightcurve at the bottom of the plot (3198 A). 
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3 HST OBSERVATIONS 

Three primary eclipses of SV Cam were observed by the 
HST, using the Space Telescope Imaging Spectrograph with 
the G430L grating. The observations used 9 spacecraft orbits 
and spanned 5 days at 2 day intervals from 1-5 November 
2001 as shown in Table 1. Summing the recorded counts 
over the observed wavelength range 2900 A to 5700 A yields 
a photometric lightcurve. Figure^ shows the photometry of 
the 3 eclipses observed during the 9 orbits. The observations 
have a cadence of 40s and a photometric precision of 0.0002 
magnitudes (S:N 5000) per 30 s exposure. The observations 
and t he data reduction m ethod are explained in greater de- 
tail in I Jeffers etafl (120051) . 



4 WAVELENGTH DEPENDENCE OF LIMB 
DARKENING 

The variation of the HST lightcurve with wavelength was 
determined by phasing the three primary eclipses together 
and dividing each spectrum into 10 bands of equal flux, as 
shown in Figure [3] The increased fluctuations at shorter 
wavelengths are consistent with the presence of signatures 
bright magnetic activity, as they are strongest in the bluest 
wavelength band. The fluctuations at blue wavelengths are 
not a short-timescale phenomenon, but rather a mismatch 
in flux levels arising from a change in the star's UV flux 
from second to third visit. 

The curvature of the eclipse profile between second 
and third contact illustrates clearly that limb darkening in- 
creases to wards shorter wavelengths. The wavelength de- 
pendence ofljmb-djxkening is revealed by the models col- 
lected bv IClaretl (120001) . but more fundamentally by direct 
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Model Grid 


Log g 


L/H 


Fitted 


Spotted 


Primary 


Secondary 


Polar Spot 








Temp (K) 


Temp. (K) 


Radius (tq) 


Radius (rQ) 


(degrees) 


PHOENIX 


4.5 


2.0 


6038±58 


5935 ±28 K 


1.235 ± 0.003 


0.727 ± 0.003 


46.5 ± 8 


ATLAS 


4.5 


1.25 


5972±59 


5840 ±53 K 


1.241 ± 0.003 


0.729 ± 0.002 


45.7 ± 9 



Table 2. Geometric binary system parameters computed for SV Cam with a reduced photospheric temperature due to the star being 
peppered with small spots, and increased radii due to the presence of a polar cap. L/H is the mixing-length to pressure scale height ratio 
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Figure 4. HST lightcurves with the longest (5596 A) and shortest 
wavelengths (3198 A), showing the variation of limb darkening 
with wavelength and the temporal evolution of bright magnetic 
features over the 5 day timespan our observations. The variation 
in the 3198A band is of the order of 5%. 

observations of the Sun, for example iNeckel fc Labsl il994ll . 
Both the brightness variation and the variation of the limb 
darkening with wavelength are illustrated by plotting the 
bluest and reddest wavelength bands, as shown in Figure 0] 



5 FIRST DERIVATIVE PROFILES OF 
LIGHTCURVES 

5.1 Basic properties 

As the primary star dominates the light of the binary sys- 
tem, the cooler secondary acts as a dark occulting disc 
scanning across the equatorial region of the primary during 
eclipse. The variation of the specific intensity as a function 
of limb angle across the primary star shows the degree of 
limb darkening. 

To determine the optimally fitting model atmosphere to 
the HST lightcurves it is necessary to examine the degree of 
curvature in the eclipse profile. We use the first derivative, 
with respect to phase, of the lightcurves to determine the 
rate of change of the eclipsed flux with phase. This is the 
first time that this method has been applied to photometric 
data of an eclipsing binary system. The numerical derivative 
is computed using the Interactive Data Language (IDL) rou- 
tines deriv and derivsig which employ 3-point Lagrangian 
interpolation. 

The profiles of the first derivative for each of the 10 HST 
wavebands are shown in Figure The first contact point is 
at phase 0.906±0.003, the second contact point is at phase 
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Figure 5. First derivative profiles for each of the 10 HST 
lightcurves as plotted in Figure [3] The presence of a bright mag- 
netic feature is the cause of the variation in waveband 3198A. 



0.977±0.003, the third contact point is at 1.023±0.003, and 
the fourth contact point is at 1.093±0.003. The gradient of 
the first derivative profile between second and third con- 
tact points is indicative of the degree of limb darkening in 
the lightcurve. In the first derivative of the model with no 
limb darkening, this region the profile is flat (FigureHjJ. The 
variation of limb darkening with wavelength is illustrated 
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Figure 6. First derivative profiles, constructed with the same 
stellar parameters as in Figure |^] but with no limb darkening. 
The slight negative slope during total eclipse results from a non- 
spherical star caused by tidal interactions in the binary system. 



Figure 8. First derivative profiles for the longest (5596 A) and 
shortest wavelengths (3198 A), as plotted in Figure^] using inten- 
sities interpolated from the PHOENIX grid at the best fit effective 
temperature and gravity shown in Table 2. 
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Figure 7. First derivative profiles of the observed lightcurve of 
wavelengths 5596 A and 3707 A for total eclipse, i.e. between sec- 
ond and third contact, showing the variation of limb darkening 
with wavelength. 



in Figure [7| where the first derivative lightcurves from the 
longest wavelength (5596 A) and the second shortest wave- 
length (3707 A) are plotted between the second and third 
contact points. The first derivative profile of the shortest 
wavelength contains too much deviation, caused by a bright 
feature on the primary's surface, to provide a clear example. 

The variation of the first derivative profile with wave- 
length is also visible in the profile using the phoenix model 
atmosphere, as shown in Figure |H| In contrast to this Fig- 
ure|U]shows the same first derivative profiles based on the AT- 
LAS model atmosphere. In the ATLAS first derivative profiles 
there is more variation between the two wavelength bands 
than in the phoenix model atmosphere. The variation with 
wavelength relates directly to the variation of specific inten- 
sity with wavelength as previously shown in Figure 
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Figure 9. First derivative profiles for with the longest (5596 A) 
and shortest wavelengths (3198 A), using intensities interpolated 
from the ATLAS grid at the best fit effective temperature and 
gravity shown in Table 2. 



5.2 First derivative lightcurve fitting 

The first derivative profiles of the HST observations, in 
each of the 10 wavelength bands (Figure |5J , are fitted us- 
ing reduced \ 2 comparison fits to the first derivatives of the 
phoenix and ATLAS model atmospheres of the same wave- 
length as the observations. In this analysis the larger error 
bars are those of the models, given the large uncertainties 
shown in Table 2, and the finite number of planer elements 
used in the modelling of the lightcurves. We estimate these 
errors to be of the order of 1%. 

The fits to the observations are determined using (i) 
the totally eclipsed section of the lightcurve between sec- 
ond and third contact points only, and (ii) the region of 
the eclipse between the first and fourth contact. Figure [TCI 
and Figure ITT1 respectively show the results for case (i) and 
(ii) at the wavelength 4535 A. The \ 2 values for all of the 
lightcurves are summarised in Table [3] 

The results show that between first and fourth con- 
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Wavelength phoenix(I) ATLAS(l) phoenix(2) ATLAS (2) 



3198 


367.27 


374.69 


285.64 


286.85 


3707 


8.41 


8.82 


7.73 


7.66 


4052 


5.82 


7.19 


6.12 


6.39 


4313 


5.70 


5.49 


4.90 


5.26 


4539 


4.84 


4.93 


4.81 


5.11 


4753 


5.00 


4.21 


5.04 


4.81 


4952 


6.08 


5.57 


4.56 


4.83 


5169 


4.60 


5.56 


4.53 


5.11 


5384 


4.50 


4.55 


4.91 


5.25 


5598 


5.09 


5.53 


6.62 


7.06 



Table 3. Best fitting reduced x 2 values for lightcurves fitted in 
the regions in the primary eclipse profile between (1) first and 
fourth contact, and (2) second and third contact. This table is 
graphically shown in Figures \l'A\ 
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Figure 11. The first derivative profile centred at 4539A, plotted 
with the PHOENIX and ATLAS model atmospheres. The phase range 
is between first and fourth contact. 
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Figure 10. The first derivative profile centred at 4539A, plotted 
with the PHOENIX and ATLAS model atmospheres. The phase range 
is between second and third contact. The lower plot shows the 
residual (observations - model) of this fit. 



tact phoenix gives the best fit, except for the wavelength 
band centred at 3198A and 3707A. However, the best fitting 
model in the region of total eclipse between second and third 
contact shows a wavelength dependence. For wavebands cen- 
tred at 4313A, 4539A, and 4743A ATLAS provides the best 
fit, while phoenix is the best fitting model atmosphere at 
shorter and longer wavelengths. 
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Figure 12. The second derivative profile using the complete HST 
lightcurve, plotted with the PHOENIX and ATLAS model atmo- 
spheres. 



6 SECOND DERIVATIVE LIGHTCURVE 
FITTING 

As with taking the first derivative of the lightcurve to de- 
termine the rate of change of the slope in the eclipse profile, 
we now determine the rate of change of the first derivative 
profile, i.e. the second derivative of the eclipse profile. The 
signal-to-noise ratio of the observed data set was insufficient 
to determine useful second derivative profiles for each of the 
10 HST lightcurves. Instead we take the second derivative 
of the complete HST lightcurve and fit phoenix and AT- 
LAS model atmospheres centred at 4670 A. We fit the second 
derivative profile between first and fourth contact as between 
second and third contact we would only fit numerical noise. 

The best fitting model atmosphere is phoenix with a 
relative \ 2 of 7.02, while the ATLAS model atmosphere has a 
relative \ 2 of 8.73. The fit to the second derivative profile for 
the phoenix and the ATLAS model atmospheres are shown 
in Figure WI\ The jitter in the models is caused by numer- 
ical noise arising from the finite number of planar surface 
elements used to model the star in the synthesis code. 
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Figure 13. The variation of reduced \ 2 as a function of wave- 
length for phases between second and third contact (top) and first 
and fourth contact (bottom) as tabulated in Table l3l 



7 DISCUSSION AND CONCLUSIONS 

The best-fitting geometric parameters determined using 
phoenix and ATLAS model atmospheres are in good agree- 
ment. As shown in Figure Q the predicted specific intensity 
for the ATLAS models is greater at the limb than for phoenix 
models, which will consequently make the primary star big- 
ger. To compensate for this the best-fitting temperature of 
the primary star is slightly cooler than fitted using phoenix 
models. In this analysis we solved the geometric system pa- 
rameters of the lightcurve using complete lightcurve rather 
than solving the parameters individually for each of the 10 
sub-lightcurves. From Figure we would expect the best- 
fitting radii, solved using phoenix and ATLAS, to be closer 
at redder wavelengths than at bluer wavelengths, but would 
not differ enough to alter the conclusions of this analysis. 

We have shown the wavelength dependence of limb 
darkening by sub-dividing the HST lightcurve into 10 bands 
of equal flux. The variation of flux between first and fourth 
contact shows that the limb darkening decreases towards 
longer wavelengths, co nfirmin g publ ished limb darkening 
values, for example by IClareji (120001) . and as observed on 
the Sun Necke]^^L i abs]j^9J) and in terferometrically in K- 
giants by iMozurkewich et al.l (120031) . The splitting of the 
HST lightcurve into 10 wavelength bands also highlights the 
presence of a time variable bright bright feature, possibly an 
active region or plage on the surface of the primary star, vis- 
ible in the bluest wavelength band (3198A). The temporal 
variation of the bright feature is of the order of 3 days as 
it comes into view on the primary star between the second 
and third HST visit. 

The ratio of the temperatures of the two stars has the 
effect that the secondary star acts as a dark occulting disc 
that scans the surface of the primary star. During partial 
eclipse (i.e. between first and second, and third and fourth 
contacts) the curvature of the lightcurve provides informa- 



tion about the variation of the specific intensity with limb 
angle of the primary star. During total eclipse the curva- 
ture of the lightcurve provides information on the variation 
of specific intensity as a function of limb angle. The first 
derivative profile for each of the 10 HST wavelength bands 
clearly indicates the change in slope as a function of phase. 
Figure [7| shows the wavelength variation of the gradient of 
the first derivative profile of the HST lightcurves centred at 
3707A and 5596A. The slope of the shorter wavelength is 
steeper than for the longer wavelength, consistent with the 
limb darkening decreasing towards longer wavelengths. 

The best fitting model atmosphere is determined using 
X 2 comparison fit. The first derivative profile of the mod- 
elled lightcurves, generated with phoenix and ATLAS model 
atmospheres, was fitted to the HST lightcurves. The results 
show that the majority of the results differ by less than la 
making the differences between the two models largely in- 
significant. 

Surface brightness reconstruction techniques such as 
Doppler imaging, and eclipse mapping rely on the informa- 
tion content of surface areas wit h different distan c es fro m 
the rotation axis (see review bv ICollier Cameron! (|2001)). 
To detect starspots Doppler imaging uses the relative in- 
tensity contributions, calculated from model atmospheres, 
to represent the different surface elements. To reconstruct 
an accurate surface brightness distribution it is essential 
to know how parameters, such as the limb darkening, can 
alter the intensity values across the stellar disc. To date 
there are many surface brightness images reconstructed us- 
ing Doppler imaging and eclipse mapping techniques on 
stars with similar s pectral types to SV Cam. Examples in- 
clude: He699 G2V jjeffers et al.ll2002l). AE Phe G0V+ F8V 
jBarnes et alJl2004l). Lu Lup, G 2V JPonati et ai]l200fj) and 
R58. G2V iMarsdaaet al .120051) . reconstructed using ATLAS 
plane-parallel model atmospheres. The lightcurve modelling 
results of this work clearly show that there is no distin- 
guishable difference between the two models using the high 
signal-to-noise SV Cam observations. 

The first derivative profiles show a small excess in the 
observed flux at phase w 0.9825 (just after second contact), 
compared with the fitted model atmospheres. In contrast, 
there is a slight decrease in the observed flux at phase 1.015, 
i.e. just before third contact. The increase and decrease of 
the first derivative at these points could be evidence for an 
additional emission. As this light excess is located just be- 
fore the second contact point, and the reverse just before 
the third contact point, it could indicate that the very edge 
of the secondary's limb is transparent to the light of the 
primary star. 

The signal-to-noise ratio of this data set was not high 
enough to determine useful second derivative profiles for 
each of the 10 HST lightcurves. Instead we take the second 
derivative of the complete HST lightcurve and fit phoenix 
and ATLAS model atmospheres centred at 4670A to the ob- 
served lightcurve. The jitter in the models is caused by nu- 
merical noise arising from the finite number of planar surface 
elements used to model the star in the synthesis code. We fit 
the lightcurve between first and fourth contact as there is in- 
sufficient structure to fit between second and third contact. 
The results show that the phoenix model atmosphere code 
gives a marginally better fit at the limb of the star. However, 
phoenix does not provide an exact fit, which could indicate 
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that the observed cut-off in the lim b intensity is steepe r than 
predicted. This could explain whv ljeffers et all <200Sh could 
not completely remove the strong discontinuities in the ob- 
served minus computed residual. 
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APPENDIX A: DETERMINATION OF THE 
PERCENTAGE OF SPOT COVERAGE AND 
POLAR CAP TEMPERATURE FOR ATLAS 
MODEL ATMOSPHERES 

In this appendix we describe the method used to determine 
the reduced photospheric temperature due to the presence 
of many unresolvable spots and the polar cap size, as quoted 
in Tab le Bl for the ATLAS model atmospheres. In a related 
paper, Ijeffers et al.l <|2005T) showed that in order to fit the 
HST lightcurve of SV Cam it was necessary to reduce the 
photospheric temperature, to mimic the peppering of the 
primary star's surface with small starspots, and to include 
a polar cap. In this appendix we determine the best-fitting 
lightcurve solution to the SV Cam lightcurve using ATLAS 
model atmosphere. It is important to do this to not to in- 
troduce an inherent bias to the results of this paper. 

Al Unresolved spot coverage 

In this section deter mine the unresolved spot coverage, fol- 
lowing the method of I Jeffers et alJ a2005f) . We assume that 
the unresolved spot coverage comprises many small spots 
peppering the primary star, and a polar cap which is not 
possible to resolve from a photometric lightcurve. 

A 1.1 Temperature Fitting 

In the method of lJeffers et alJ j2005t) the combination of the 
SV Cam lightcurve and the Hipparcos parallax is used to de- 
termine the primary and secondary temperatures. Knowing 
the radii of the two stars we can evaluate the flux contribu- 
tion from the secondary star relative to that of the primary 
star. The best-fitting combination of primary and secondary 
temperatures is determined using \ 2 minimisation, where a 
scaling factor 7 is included to ensure that the shape of the 
spectrum is fitted rather than the absolute flux levels. The 
resulting x 2 landscape plot is shown in Figure lA"Tl The min- 
imum value occurs at 5973±31 K and 4831±103K for the 
primary and secondary stars respectively. 

The temperature of the primary star is then determined 
by isolating the primary star's spectrum. To achieve this we 
subtracted a spectrum outside of the primary eclipse from 
one inside of the eclipse which results in the spectrum of 
the primary star but with the radius of the secondary star. 
The temperature is fitted using \ 2 minimisation, where the 
minimum temperature is determined by a parabolic fit. This 
results in a minimum primary temperature of 5872±59 K 
(Figure fA2ll . with the errors determined by setting Ax 2 =l- 

A 1.2 Fractional Starspot Coverage 

Following the conclusions o flJeffers" et al.H2005D we attribute 
the missing flux to be indicative of small unresolvable spots 
on the primary star's surface. The dark starspot filling factor 
is given by: 

a = 1-7 (Al) 

where a is the fractional starspot coverage and 7 is 
the scaling factor. The interpolated scaling factor for the 
primary temperature is determined as shown in Figure lA"3l 
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Primary Radius 



Figure Al. Contour plot of the x landscape of the combined 
primary and secondary stars. The minimum x 2 value occurs at 
5973 ± 31 K, and 4831 ± 103 K, for the primary and secondary 
stars respectively. From the centre of the plot the first contour 
ellipse represents the 1 parameter 1 <r confidence limit at 63.8%, 
the second ellipse represents the 2 parameter 1 a confidence limit 
at 63.8% whilst the third ellipse represents the 2 parameter 2.6 a 
99% confidence limit. 
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Figure A2. Parabolic fit to determine the primary temperature 
minimum to be 5872±53 K 



The scaling factor is 0.76, resulting in a fractional coverage 
of dark starspots of 24%. 



A 1.3 Polar Cap 

The determination of the spot coverage fraction only ac- 
counts for the flux deficit in the eclipsed equatorial lati- 
tudes of the primary star. Extending the 24% spot cover- 
age to_the_ent^ surface of the primary star (as described 
by j.Ieffers et al.l l2"005T)) we find that there is an additional 
13.5% flux deficit. The binary eclipse-mapping code DoTs 
is used to model artificial polar spots on the surface of SV 
Cam, to include effects resulting from the star being a sphere 
and not a disc, limb and gravity darkening and spherical 
oblateness. We model the fractional decrease in stellar flux 
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Figure A3. Scaling factor (7) as a function of Primary Temper- 
ature (atlas models). The scaling factor for T=5872K is deter- 
mined by a quadratic fit to these points. 



DoTs Model Flux Decrease 
SV Cam (ATLAS) 
SV Cam (upper limit) 
SV Cam (lower limit) 
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Figure A4. Fractional decrease in the stellar flux of SV Cam as 
a function of theoretical polar cap size 



as a function of polar spot size and determine the polar spot 
radius to be 43.5±6° (Fig. lA4)l . 



A2 Lightcurve Modelling 

To compare ATLAS and PHOENIX model atmospheres we 
need to determine the best-fitting binary system parameters 
to the observed SV Cam lightcurve for each model atmo- 
sphere separately to avoid introducing a an inherent bias 
to our results. We include the presence of high unresolvable 
spot coverage and polar caps in t he lightcurve fit by using 



50 



cap s m t n 

the method of lJeffers et alJ 12005). Such spot coverage has 
been shown by that method to have a significant impact on 
the binary system parameters. 

A2.1 Reduced Photospheric Temperature 

The peppering of small starspots poses a limitation on im- 
age reconstruction techniques from photomet ric lightcurve s 
such as Maximum Entropy eclipse mapping jJeffersll2005Tl . 
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Figure A5. The decrease of the primary star's apparent pho- 
tospheric temperature as a function of the percentage of spot 
coverage on its surface. 
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Figure A6. Quadratic fit to the variation of x 2 as a function of 
polar spot size. 



To model the presence of many dark unresolvable starspots 
we decrease the apparent photospheric temperature of the 
star. To determine the reduction in the photospheric temper- 
ature of the star we model starspot distributions equating 
to 1.8%, 6.1%, 18%, 48% and 100% of the stellar surface on 
an immaculate SV Cam. For each model the initial photo- 
spheric temperature is 5904 K and the spot temperature is 
4400 K. Each of these starspot distributions is modelled as a 
photometric lightcurve. Using the Maximum Entropy eclipse 
mapping method we determine the best-fitting temperature 
to each model lightcurve using a \ 2 grid-search method. A 
quadratic fit to the best-fitting temperatures show that for 
a starspot coverage of 24%, the apparent photospheric tem- 
perature is 5840 K f Figure IXSY 

A2.2 Polar Cap 

We include a polar cap in our analysis to verify the results in 
the pr evious section, following the method of I.Teffers et alJ 
(2005). The polar spot is assumed to be at 4500 K, circu- 
lar, centred at the pole, and is in addition to the peppered 
spot distribution as described above. For each polar spot 
size (from 40° to 50°) the minimum primary and secondary 
radii are determined using a \ 2 contour map. These mini- 
mum \ 2 values are plotted as a function of polar spot size 
in Figure lA"6l The best fitting polar cap size, 45.7° was de- 
termined from the minimum of a quadratic function fitted 
to these points. The grid search of radii is repeated using 
a fixed value for the polar spot size. The results for the \ 2 
minimisation are summarised in Table [2] shown as a contour 
map in Figure lATI 




1.235 1.24 1.245 

Primary Radius 



Figure A7. Contour plot of the \ 2 landscape for the primary 
and secondary radii with a 45.7° polar spot. From the centre of 
the plot, the first contour ellipse is the 1 parameter lrj confidence 
limit at 63.8%, the second contour ellipse is the 2 parameter la 
confidence limit at 63.8%, and the third contour ellipse is the 2 
parameter 2.6c confidence limit at 99%. 

phoenix and ATLAS model atmospheres at the F9V+K4V 
spectral types. The results are summarised in Table [3] 

This paper has been typeset from a TjEX/ I^TeX file prepared 
by the author. 



A3 Summary 

We have shown in this section that the best fitting binary 
system parameters determined using the ATLAS model at- 
mosphere are in agreement with those determ ined using 
phoenix model atmospheres (fjeflers et alj l2005l. This fur- 
ther shows that at there is no significant difference between 
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ABSTRACT 

The variation of the specific intensity across the stellar disc is essential input 
parameter in surface brightness reconstruction techniques such as Doppler imaging, 
where the relative intensity contributions of different surface elements are important 
in detecting starspots. We use phoenix and ATLAS model atmospheres to model 
lightcurves derived from high precision (S/N ~ 5000) HST data of the eclipsing binary 
SV Cam (F9V + K4V), where the variation of specific intensity across the stellar disc 
will determine the contact points of the binary system lightcurve. For the first time 
we use x 2 comparison fits to the first derivative profiles to determine the best-fitting 
model atmosphere. We show the wavelength dependence of the limb darkening and 
that the first derivative profile is sensitive to the limb-darkening profile very close 
to the limb of the primary star. It is concluded that there is only a marginal differ- 
ence (< lex) between the x 2 comparison fits of the two model atmospheres to the HST 
lightcurve at all wavelengths. The usefulness of the second derivative of the light-curve 
for measuring the sharpness of the primary's limb is investigated, but we find that the 
data are too noisy to permit a quantitative analysis. 

Key words: 

stars: activity, stars: spots binaries: eclipsing, stars:atmospheres,methods:numerical 



1 INTRODUCTION 

Limb darkening effects in stellar atmospheres have impor- 
tant implications throughout stellar astrophysics where a de- 
termination of the surface brightness distribution is impor- 
tant. Recent work using Doppler imaging and micro- lensing 
events has shown that commonly-used analytical limb dark- 
ening laws fail to match stellar observations at the limb of 
the star (??). Other micro-lensing results (?) show that the 
intensity predictions from model atmospheres are discrepant 
in the case of a K-giant. 

Surface brightness reconstruction techniques such as 
Doppler imaging and eclipse mapping rely on the informa- 
tion content of surface areas with different distances from 
the rotation axis (see review by ?). To detect starspots 
Doppler imaging uses the relative intensity contributions, 
calculated from model atmospheres, of the different surface 
elements. To reconstruct an accurate surface brightness dis- 
tribution it is essential to know how parameters, such as 



the limb darkening, can alter the intensity values across the 
stellar disc. 

High inclination eclipsing binary systems can be used as 
probes to determine the variation of specific intensity at the 
stellar limb. If limb darkening showed a smooth transition in 
specific intensity at the limb of the star, the contact points 
of eclipses would appear less abrupt and slightly displaced in 
phase relative to models with limb darkening laws derived 
from plane-parallel atmospheres, where the cutoff is very 
sharp. The sharp cutoff in plane parallel atmospheres results 
from the optical depth of the rays being infinite at the limb. 

In November 2001 we were awarded 9 orbits of 
HST/STIS time to eclipse-map the inner face of the F9V 
primary of the totally eclipsing binary SV Cam. SV Cam 
(F9V + K4V) is a synchronously rotating RS CVn binary 
with a period of 0.59d. We obtained spectrophotometric 
lightcurves of 3 primary eclipses with a signal-to-noise ratio 
of 5000. The first analysis of these data, by ?, determined the 
radii of the primary and secondary stars. When the resulting 
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lightcurve was subtracted from the observed data, the resid- 
ual lightcurve showed strong peaks at phases of contact. ? 
then showed that these mismatches are reduced significantly, 
but not eliminated, when a polar cap and a reduction in the 
photospheric temperature, to synthesise high spot coverage, 
are imposed on the image. 

As there is a significant temperature difference between 
the primary and secondary stars, the secondary star acts as 
a dark occulting disc as it eclipses the primary star. The 
variation of brightness as a function of phase reflects the 
degree of limb darkening on the primary star. In this paper 
we determine the best fitting model atmosphere by fitting 
the models to the brightness variations as the secondary 
scans the inner face of the primary star, using the first and 
second derivatives of the HST in 10 wavelength bands. We 
discuss the implications these results have on results from 
Doppler Imaging. 



2 MODEL ATMOSPHERES 

In this paper, two well established stellar atmosphere codes 
are used; the phoenix model atmosphere code (?) that 
uses spherical atmospheres and the ATLAS model atmosphere 
code (?) that uses plane parallel atmospheres. 

2.1 ATLAS 

We use the plane-parallel ATLAS9 model atmospheres from 
the Kurucz CD-ROMS (?). We integrate the intensity val- 
ues over the wavelength range of our observations 2900A to 
5700A. We use temperature models from 3500 K to 6500 K, 
with 250 K interval, across 17 limb angles. The limb angle [i 
is defined by fi— cos 9, where 9 is the angle between the line 
of sight and the normal vector of the local surface element. 
The treatment of convection is based on the mixing length 
theory with approximate overshoot (ref. ?) with a mixing 
length to scale height ratio of 1.25. The variation of specific 
intensity, i.e. at ^i=l,as a function of wavelength and limb 
angle is shown in Figure 

2.2 PHOENIX 

The general input physics set-up of the phoenix model at- 
mosphere code is discussed in ?. The main advantage of 
using this code is that it is based on spherical geometry 
(spherical radiative transfer) LTE rather than traditional 
plane-parallel structure. NLTE effects are considered to be 
insignificant in this application. 

The synthetic spectra are based on an extension of 
the grid of PHOENIX model atmospheres described by ?. 
This extended grid includes surface gravities larger than 
log(<?) > 3.0 needed for main sequence stars. These models 
are as described by ?, but include an updated molecular line 
list. The models are computed in spherical geometry with 
full atomic and molecular line blanketing using solar elemen- 
tal abundances. In these models, the stellar mass is 0.5 Mq 
and the convection treatment assumes a mixing-length to 
pressure scale height ratio of 2 with no overshooting. There 
are 117 synthetic spectra in total. The effective temperature 
runs from 2700K to 6500 K in 100K steps at three surface 
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Figure 1. The predicted variation in specific intensity across the 
stellar disc as a function of wavelength for the longest (5596A) 
and the shortest (3198A) wavelengths in this analysis (T=6000K, 
log g=4.5). 

gravities: log(g) = 4.0, 4.5, and 5.0. The wavelength resolu- 
tion of these synthetic spectra is lA. 

The variation of specific intensity as a function of limb 
angle is shown in Figure Q The difference between the 
phoenix and the ATLAS model atmospheres at the limb of 
the star results primarily from the effects of spherical geom- 
etry of the phoenix model. In models using spherical geom- 
etry, there is a finite optical depth for rays close to the limb, 
while with plane parallel models the optical depth of the rays 
is infinite at the limb providing significant intensities down 
to /J.—0. Other contributing effects include overshooting and 
the mixing length ratio. 

In addition to the limb effects, Fig.^shows that ATLAS 
intensity profiles are overall brighter than phoenix intensity 
profiles. We attribute this difference to the use of overshoot- 
ing in the ATLAS models. ATLAS models with overshooting 
have been shown to have weaker limb darkening in the blue 
relative to models without overshooting, for instance, over- 
shooting models have been shown to better fit solar limb 
darkening observations than non-overshooting models ?. A 
similar result holds for Procyon (F5 IV- V); ref. ? for a de- 
tailed study of the effects of 1-D and 3-D convection on limb 
darkening. 

2.3 Model lightcurves 

We use the eclipse-mapping code DoTS (?) to model the pri- 
mary eclipse lightcurve for each model atmosphere, and to 
determine which model atmosphere best fits the complete 
HST lightcurve. The input data to each model comprises: 
the variation of the specific intensity as a function of limb 
angle, as shown in FigureQfor both model atmospheres con- 
sidered; the primary and secondary radii; a reduced primary 
photospheric temperature and polar spot. Gravity darkening 
is also included according the description of ?. 

We include a reduced photospheric temperature to ac- 
count for a stellar surface that is peppered with small spots 
which are too small to be resolved through eclipse mapping. 
Following the results of ? we include a polar cap in the mod- 
elled lightcurve to optimise the fit of the lightcurve to the 
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Visit 


Obs. Date 


UT 


UT 


Exposure 


No of 






Start 


End 


Time(s) 


Frames 


1 


01 November 2001 


20:55:56 


01:00:17 


30 


165 


2 


03 November 2001 


14:34:29 


18:38:21 


30 


165 


3 


05 November 2001 


09:49:17 


13:52:22 


30 


165 



Table 1. HST Observations of SV Cam, where the exposure time 
is per frame. 
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Figure 2. The comprising sections of the primary eclipse from 
each of the 3 HST visits. The flux density has been summed 
between 2900 A and 5700 A. An off-set of 4xl0 6 ergs/s/ cm has 
been added for clarity. 



data. The method for modelling the photometric lightcurve 
to include a reduced photospheric temperature and a polar 
cap is described in Appendix A for the ATLAS model atmo- 
sphere. The binary system parameters are summarised in 
Table |3 The lightcurve solutions for the two model atmo- 
spheres are in good agreement. 



3 HST OBSERVATIONS 

Three primary eclipses of SV Cam were observed by the 
HST, using the Space Telescope Imaging Spectrograph with 
the G430L grating. The observations used 9 spacecraft orbits 
and spanned 5 days at 2 day intervals from 1-5 November 
2001 as shown in Table 1. Summing the recorded counts 
over the observed wavelength range 2900 A to 5700 A yields 
a photometric lightcurve. Figure |21 shows the photometry of 
the 3 eclipses observed during the 9 orbits. The observations 
have a cadence of 40s and a photometric precision of 0.0002 
magnitudes (S:N 5000) per 30 s exposure. The observations 
and the data reduction method are explained in greater de- 
tail in ?. 



4 WAVELENGTH DEPENDENCE OF LIMB 
DARKENING 

The variation of the HST lightcurve with wavelength was 
determined by phasing the three primary eclipses together 
and dividing each spectrum into 10 bands of equal flux, as 
shown in Figure [3] The increased fluctuations at shorter 
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Figure 3. The variation of the HST lightcurve with wavelength. 
Each lightcurve has been plotted with an offset of 5 X 10~ 12 from 
the lightcurve at the bottom of the plot (3198 A). 



wavelengths are consistent with the presence of signatures 
bright magnetic activity, as they are strongest in the bluest 
wavelength band. The fluctuations at blue wavelengths are 
not a short-timescale phenomenon, but rather a mismatch 
in flux levels arising from a change in the star's UV flux 
from second to third visit. 

The curvature of the eclipse profile between second 
and third contact illustrates clearly that limb darkening in- 
creases to wards shorter wavelengths. The wavelength de- 
pendence of limb-darkening is revealed by the models col- 
lected by ?, but more fundamentally by direct observations 
of the Sun, for example ?. Both the brightness variation and 
the variation of the limb darkening with wavelength are illus- 
trated by plotting the bluest and reddest wavelength bands, 
as shown in Figure 0] 
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Model Grid 


Log g 


L/H 


Fitted 


Spotted 


Primary 


Secondary 


Polar Spot 








Temp (K) 


Temp. (K) 


Radius (tq) 


Radius (rQ) 


(degrees) 


PHOENIX 


4.5 


2.0 


6038±58 


5935 ±28 K 


1.235 ± 0.003 


0.727 ± 0.003 


46.5 ± 8 


ATLAS 


4.5 


1.25 


5972±59 


5840 ±53 K 


1.241 ± 0.003 


0.729 ± 0.002 


45.7 ± 9 



Table 2. Geometric binary system parameters computed for SV Cam with a reduced photospheric temperature due to the star being 
peppered with small spots, and increased radii due to the presence of a polar cap. L/H is the mixing-length to pressure scale height ratio 
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Figure 4. HST lightcurves with the longest (5596 A) and shortest 
wavelengths (3198 A), showing the variation of limb darkening 
with wavelength and the temporal evolution of bright magnetic 
features over the 5 day timespan our observations. The variation 
in the 3198A band is of the order of 5%. 

5 FIRST DERIVATIVE PROFILES OF 
LIGHTCURVES 

5.1 Basic properties 

As the primary star dominates the light of the binary sys- 
tem, the cooler secondary acts as a dark occulting disc 
scanning across the equatorial region of the primary during 
eclipse. The variation of the specific intensity as a function 
of limb angle across the primary star shows the degree of 
limb darkening. 

To determine the optimally fitting model atmosphere to 
the HST lightcurves it is necessary to examine the degree of 
curvature in the eclipse profile. We use the first derivative, 
with respect to phase, of the lightcurves to determine the 
rate of change of the eclipsed flux with phase. This is the 
first time that this method has been applied to photometric 
data of an eclipsing binary system. The numerical derivative 
is computed using the Interactive Data Language (IDL) rou- 
tines deriv and derivsig which employ 3-point Lagrangian 
interpolation. 

The profiles of the first derivative for each of the 10 HST 
wavebands are shown in Figure [S] The first contact point is 
at phase 0.906±0.003, the second contact point is at phase 
0.977±0.003, the third contact point is at 1.023±0.003, and 
the fourth contact point is at 1.093±0.003. The gradient of 
the first derivative profile between second and third con- 
tact points is indicative of the degree of limb darkening in 
the lightcurve. In the first derivative of the model with no 
limb darkening, this region the profile is flat (Figure The 
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Figure 5. First derivative profiles for each of the 10 HST 
lightcurves as plotted in Figure [3] The presence of a bright mag- 
netic feature is the cause of the variation in waveband 3198A. 



variation of limb darkening with wavelength is illustrated 
in Figure |7] where the first derivative lightcurves from the 
longest wavelength (5596 A) and the second shortest wave- 
length (3707 A) are plotted between the second and third 
contact points. The first derivative profile of the shortest 
wavelength contains too much deviation, caused by a bright 
feature on the primary's surface, to provide a clear example. 
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Figure 6. First derivative profiles, constructed with the same 
stellar parameters as in Figure |^] but with no limb darkening. 
The slight negative slope during total eclipse results from a non- 
spherical star caused by tidal interactions in the binary system. 



Figure 8. First derivative profiles for the longest (5596 A) and 
shortest wavelengths (3198 A), as plotted in Figure^] using inten- 
sities interpolated from the PHOENIX grid at the best fit effective 
temperature and gravity shown in Table 2. 
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Figure 7. First derivative profiles of the observed lightcurve of 
wavelengths 5596 A and 3707 A for total eclipse, i.e. between sec- 
ond and third contact, showing the variation of limb darkening 
with wavelength. 

The variation of the first derivative profile with wave- 
length is also visible in the profile using the phoenix model 
atmosphere, as shown in Figure [§] In contrast to this Fig- 
ure[5]shows the same first derivative profiles based on the AT- 
LAS model atmosphere. In the ATLAS first derivative profiles 
there is more variation between the two wavelength bands 
than in the phoenix model atmosphere. The variation with 
wavelength relates directly to the variation of specific inten- 
sity with wavelength as previously shown in Figure 

5.2 First derivative lightcurve fitting 

The first derivative profiles of the HST observations, in 
each of the 10 wavelength bands (Figure , are fitted us- 
ing reduced x 2 comparison fits to the first derivatives of the 
phoenix and ATLAS model atmospheres of the same wave- 
length as the observations. In this analysis the larger error 
bars are those of the models, given the large uncertainties 
shown in Table 2, and the finite number of planer elements 
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Figure 9. First derivative profiles for with the longest (5596 A) 
and shortest wavelengths (3198 A), using intensities interpolated 
from the ATLAS grid at the best fit effective temperature and 
gravity shown in Table 2. 



used in the modelling of the lightcurves. We estimate these 
errors to be of the order of 1%. 

The fits to the observations are determined using (i) 
the totally eclipsed section of the lightcurve between sec- 
ond and third contact points only, and (ii) the region of 
the eclipse between the first and fourth contact. Figure [TU1 
and Figure ITT1 respectively show the results for case (i) and 
(ii) at the wavelength 4535 A. The \ 2 values for all of the 
lightcurves are summarised in Table [3] 

The results show that between first and fourth con- 
tact phoenix gives the best fit, except for the wavelength 
band centred at 3198A and 3707A. However, the best fitting 
model in the region of total eclipse between second and third 
contact shows a wavelength dependence. For wavebands cen- 
tred at 4313A, 4539A, and 4743A atlas provides the best 
fit, while phoenix is the best fitting model atmosphere at 
shorter and longer wavelengths. 
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Wavelength phoenix(I) ATLAS(l) phoenix(2) ATLAS (2) 



3198 


367.27 


374.69 


285.64 


286.85 


3707 


8.41 


8.82 


7.73 


7.66 


4052 


5.82 


7.19 


6.12 


6.39 


4313 


5.70 


5.49 


4.90 


5.26 


4539 


4.84 


4.93 


4.81 


5.11 


4753 


5.00 


4.21 


5.04 


4.81 


4952 


6.08 


5.57 


4.56 


4.83 


5169 


4.60 


5.56 


4.53 


5.11 


5384 


4.50 


4.55 


4.91 


5.25 


5598 


5.09 


5.53 


6.62 


7.06 



Table 3. Best fitting reduced x 2 values for lightcurves fitted in 
the regions in the primary eclipse profile between (1) first and 
fourth contact, and (2) second and third contact. This table is 
graphically shown in Figures \l'A\ 




0.92 0.94 0.96 0.98 1 1.02 1.04 1.06 1.08 
Phase 



Figure 11. The first derivative profile centred at 4539A, plotted 
with the PHOENIX and ATLAS model atmospheres. The phase range 
is between first and fourth contact. 
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Figure 10. The first derivative profile centred at 4539A, plotted 
with the PHOENIX and ATLAS model atmospheres. The phase range 
is between second and third contact. The lower plot shows the 
residual (observations - model) of this fit. 



6 SECOND DERIVATIVE LIGHTCURVE 
FITTING 

As with taking the first derivative of the lightcurve to de- 
termine the rate of change of the slope in the eclipse profile, 
we now determine the rate of change of the first derivative 
profile, i.e. the second derivative of the eclipse profile. The 
signal-to-noise ratio of the observed data set was insufficient 
to determine useful second derivative profiles for each of the 
10 HST lightcurves. Instead we take the second derivative 
of the complete HST lightcurve and fit phoenix and AT- 
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Figure 12. The second derivative profile using the complete HST 
lightcurve, plotted with the PHOENIX and ATLAS model atmo- 
spheres. 



LAS model atmospheres centred at 4670 A. We fit the second 
derivative profile between first and fourth contact as between 
second and third contact we would only fit numerical noise. 

The best fitting model atmosphere is phoenix with a 
relative \ 2 of 7.02, while the ATLAS model atmosphere has a 
relative \ 2 of 8.73. The fit to the second derivative profile for 
the phoenix and the ATLAS model atmospheres are shown 
in Figure WI\ The jitter in the models is caused by numer- 
ical noise arising from the finite number of planar surface 
elements used to model the star in the synthesis code. 



7 DISCUSSION AND CONCLUSIONS 

The best-fitting geometric parameters determined using 
phoenix and ATLAS model atmospheres are in good agree- 
ment. As shown in Figure the predicted specific intensity 
for the ATLAS models is greater at the limb than for phoenix 
models, which will consequently make the primary star big- 
ger. To compensate for this the best-fitting temperature of 
the primary star is slightly cooler than fitted using phoenix 



Modelling of HST lightcurves 7 




3600 3800 4000 4200 4400 4600 4800 5000 5200 5400 5600 
Wavelength Band 



3600 3800 4000 4200 4400 4600 4800 5000 5200 5400 5600 
Wavelength Band 

Figure 13. The variation of reduced \ 2 as a function of wave- 
length for phases between second and third contact (top) and first 
and fourth contact (bottom) as tabulated in Table l3l 



models. In this analysis we solved the geometric system pa- 
rameters of the lightcurve using complete lightcurve rather 
than solving the parameters individually for each of the 10 
sub-lightcurves. From Figure we would expect the best- 
fitting radii, solved using phoenix and atlas, to be closer 
at redder wavelengths than at bluer wavelengths, but would 
not differ enough to alter the conclusions of this analysis. 

We have shown the wavelength dependence of limb 
darkening by sub-dividing the HST lightcurve into 10 bands 
of equal flux. The variation of flux between first and fourth 
contact shows that the limb darkening decreases towards 
longer wavelengths, confirming published limb darkening 
values, for example by ?, and as observed on the Sun ? 
and interferometrically in K-giants by ?. The splitting of the 
HST lightcurve into 10 wavelength bands also highlights the 
presence of a time variable bright bright feature, possibly an 
active region or plage on the surface of the primary star, vis- 
ible in the bluest wavelength band (3198A). The temporal 
variation of the bright feature is of the order of 3 days as 
it comes into view on the primary star between the second 
and third HST visit. 

The ratio of the temperatures of the two stars has the 
effect that the secondary star acts as a dark occulting disc 
that scans the surface of the primary star. During partial 
eclipse (i.e. between first and second, and third and fourth 
contacts) the curvature of the lightcurve provides informa- 
tion about the variation of the specific intensity with limb 
angle of the primary star. During total eclipse the curva- 
ture of the lightcurve provides information on the variation 
of specific intensity as a function of limb angle. The first 
derivative profile for each of the 10 HST wavelength bands 
clearly indicates the change in slope as a function of phase. 
Figure Q shows the wavelength variation of the gradient of 
the first derivative profile of the HST lightcurves centred at 
3707A and 5596A. The slope of the shorter wavelength is 



steeper than for the longer wavelength, consistent with the 
limb darkening decreasing towards longer wavelengths. 

The best fitting model atmosphere is determined using 
X 2 comparison fit. The first derivative profile of the mod- 
elled lightcurves, generated with phoenix and ATLAS model 
atmospheres, was fitted to the HST lightcurves. The results 
show that the majority of the results differ by less than la 
making the differences between the two models largely in- 
significant. 

Surface brightness reconstruction techniques such as 
Doppler imaging, and eclipse mapping rely on the informa- 
tion content of surface areas with different distances from 
the rotation axis (see review by ?). To detect starspots 
Doppler imaging uses the relative intensity contributions, 
calculated from model atmospheres, to represent the dif- 
ferent surface elements. To reconstruct an accurate surface 
brightness distribution it is essential to know how parame- 
ters, such as the limb darkening, can alter the intensity val- 
ues across the stellar disc. To date there are many surface 
brightness images reconstructed using Doppler imaging and 
eclipse mapping techniques on stars with similar spectral 
types to SV Cam. Examples include: He699, G2V (?), AE 
Phe G0V+F8V (?), Lu Lup, G2V (?) and R58, G2V (?), re- 
constructed using ATLAS plane-parallel model atmospheres. 
The lightcurve modelling results of this work clearly show 
that there is no distinguishable difference between the two 
models using the high signal-to-noise SV Cam observations. 

The first derivative profiles show a small excess in the 
observed flux at phase » 0.9825 (just after second contact), 
compared with the fitted model atmospheres. In contrast, 
there is a slight decrease in the observed flux at phase 1.015, 
i.e. just before third contact. The increase and decrease of 
the first derivative at these points could be evidence for an 
additional emission. As this light excess is located just be- 
fore the second contact point, and the reverse just before 
the third contact point, it could indicate that the very edge 
of the secondary's limb is transparent to the light of the 
primary star. 

The signal-to-noise ratio of this data set was not high 
enough to determine useful second derivative profiles for 
each of the 10 HST lightcurves. Instead we take the second 
derivative of the complete HST lightcurve and fit phoenix 
and ATLAS model atmospheres centred at 4670A to the ob- 
served lightcurve. The jitter in the models is caused by nu- 
merical noise arising from the finite number of planar surface 
elements used to model the star in the synthesis code. We fit 
the lightcurve between first and fourth contact as there is in- 
sufficient structure to fit between second and third contact. 
The results show that the phoenix model atmosphere code 
gives a marginally better fit at the limb of the star. How- 
ever, phoenix does not provide an exact fit, which could 
indicate that the observed cut-off in the limb intensity is 
steeper than predicted. This could explain why ? could not 
completely remove the strong discontinuities in the observed 
minus computed residual. 
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APPENDIX A: DETERMINATION OF THE 
PERCENTAGE OF SPOT COVERAGE AND 
POLAR CAP TEMPERATURE FOR ATLAS 
MODEL ATMOSPHERES 

In this appendix we describe the method used to determine 
the reduced photospheric temperature due to the presence 
of many unresolvable spots and the polar cap size, as quoted 
in Table for the ATLAS model atmospheres. In a related 
paper, ? showed that in order to fit the HST lightcurve of SV 
Cam it was necessary to reduce the photospheric tempera- 
ture, to mimic the peppering of the primary star's surface 
with small starspots, and to include a polar cap. In this ap- 
pendix we determine the best-fitting lightcurve solution to 
the SV Cam lightcurve using ATLAS model atmosphere. It is 
important to do this to not to introduce an inherent bias to 
the results of this paper. 



Al Unresolved spot coverage 

In this section determine the unresolved spot coverage, fol- 
lowing the method of ?. We assume that the unresolved spot 
coverage comprises many small spots peppering the primary 
star, and a polar cap which is not possible to resolve from a 
photometric lightcurve. 



A 1.1 Temperature Fitting 

In the method of ? the combination of the SV Cam 
lightcurve and the Hipparcos parallax is used to determine 
the primary and secondary temperatures. Knowing the radii 
of the two stars we can evaluate the flux contribution from 
the secondary star relative to that of the primary star. The 
best-fitting combination of primary and secondary temper- 
atures is determined using \ 2 minimisation, where a scaling 
factor 7 is included to ensure that the shape of the spectrum 
is fitted rather than the absolute flux levels. The resulting \ 2 
landscape plot is shown in Figure |A"T1 The minimum value 
occurs at 5973±31 K and 4831±103K for the primary and 
secondary stars respectively. 

The temperature of the primary star is then determined 
by isolating the primary star's spectrum. To achieve this we 
subtracted a spectrum outside of the primary eclipse from 
one inside of the eclipse which results in the spectrum of 
the primary star but with the radius of the secondary star. 
The temperature is fitted using \ 2 minimisation, where the 
minimum temperature is determined by a parabolic fit. This 
results in a minimum primary temperature of 5872±59 K 
(Figure fA2L with the errors determined by setting Ax 2 =l- 
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Figure Al. Contour plot of the \ 2 landscape of the combined 
primary and secondary stars. The minimum \ 2 value occurs at 
5973 ± 31 K, and 4831 ± 103 K, for the primary and secondary 
stars respectively. From the centre of the plot the first contour 
ellipse represents the 1 parameter 1 a confidence limit at 63.8%, 
the second ellipse represents the 2 parameter 1 a confidence limit 
at 63.8% whilst the third ellipse represents the 2 parameter 2.6 a 
99% confidence limit. 
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Figure A2. Parabolic fit to determine the primary temperature 
minimum to be 5872±53 K 



A 1.2 Fractional Starspot Coverage 

Following the conclusions of ? we attribute the missing flux 
to be indicative of small unresolvable spots on the primary 
star's surface. The dark starspot filling factor is given by: 

a = 1-7 (Al) 

where a is the fractional starspot coverage and 7 is 
the scaling factor. The interpolated scaling factor for the 
primary temperature is determined as shown in Figure lA"3l 
The scaling factor is 0.76, resulting in a fractional coverage 
of dark starspots of 24%. 
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Figure A3. Scaling factor (7) as a function of Primary Temper- 
ature (atlas models). The scaling factor for T=5872K is deter- 
mined by a quadratic fit to these points. 
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Figure A4. Fractional decrease in the stellar flux of SV Cam as 
a function of theoretical polar cap size 

A 1.3 Polar Cap 

The determination of the spot coverage fraction only ac- 
counts for the flux deficit in the eclipsed equatorial latitudes 
of the primary star. Extending the 24% spot coverage to the 
entire surface of the primary star (as described by (?)) we 
find that there is an additional 13.5% flux deficit. The bi- 
nary eclipse-mapping code DoTs is used to model artificial 
polar spots on the surface of SV Cam, to include effects re- 
sulting from the star being a sphere and not a disc, limb 
and gravity darkening and spherical oblateness. We model 
the fractional decrease in stellar flux as a function of polar 
spot size and determine the polar spot radius to be 43.5±6° 
(Fig. E4l . 

A2 Lightcurve Modelling 

To compare ATLAS and PHOENIX model atmospheres we 
need to determine the best-fitting binary system parameters 
to the observed SV Cam lightcurve for each model atmo- 



5920 r 1 

5900 - \ 1 
5880 : \ : 
5860 r X : 

5840 r \ : 

1 5820 - \. : 

Q. 

I 5800 - N \ 

5780 - + '\ : 

5760 r 

5740 ; . 

5720 ' " 

10 20 30 40 50 60 70 80 90 100 

Percentage of Spottedness 

Figure A5. The decrease of the primary star's apparent pho- 
tospheric temperature as a function of the percentage of spot 
coverage on its surface. 

sphere separately to avoid introducing a an inherent bias 
to our results. We include the presence of high unresolvable 
spot coverage and polar caps in the lightcurve fit by using 
the method of ?. Such spot coverage has been shown by that 
method to have a significant impact on the binary system 
parameters. 

A2.1 Reduced Photospheric Temperature 

The peppering of small starspots poses a limitation on im- 
age reconstruction techniques from photometric lightcurves 
such as Maximum Entropy eclipse mapping (?). To model 
the presence of many dark unresolvable starspots we de- 
crease the apparent photospheric temperature of the star. 
To determine the reduction in the photospheric tempera- 
ture of the star we model starspot distributions equating to 
1.8%, 6.1%, 18%, 48% and 100% of the stellar surface on 
an immaculate SV Cam. For each model the initial photo- 
spheric temperature is 5904 K and the spot temperature is 
4400 K. Each of these starspot distributions is modelled as a 
photometric lightcurve. Using the Maximum Entropy eclipse 
mapping method we determine the best-fitting temperature 
to each model lightcurve using a x 2 grid-search method. A 
quadratic fit to the best-fitting temperatures show that for 
a starspot coverage of 24%, the apparent photospheric tem- 
perature is 5840 K (Figure IXSY 

A2.2 Polar Cap 

We include a polar cap in our analysis to verify the results 
in the previous section, following the method of ?. The polar 
spot is assumed to be at 4500 K, circular, centred at the pole, 
and is in addition to the peppered spot distribution as de- 
scribed above. For each polar spot size (from 40° to 50°) the 
minimum primary and secondary radii are determined using 
a x 2 contour map. These minimum \ 2 values are plotted as 
a function of polar spot size in Figure IA6I The best fitting 
polar cap size, 45.7° was determined from the minimum of 
a quadratic function fitted to these points. The grid search 
of radii is repeated using a fixed value for the polar spot 
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Figure A6. Quadratic fit to the variation of x 2 as a function of 
polar spot size. 
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Figure A7. Contour plot of the \ 2 landscape for the primary 
and secondary radii with a 45.7° polar spot. From the centre of 
the plot, the first contour ellipse is the 1 parameter la confidence 
limit at 63.8%, the second contour ellipse is the 2 parameter lrr 
confidence limit at 63.8%, and the third contour ellipse is the 2 
parameter 2.6c confidence limit at 99%. 



size. The results for the \ 2 minimisation are summarised in 
Table |3 shown as a contour map in Figure 1X71 



A3 Summary 

We have shown in this section that the best fitting binary 
system parameters determined using the ATLAS model at- 
mosphere are in agreement with those determined using 
phoenix model atmospheres (?). This further shows that 
at there is no significant difference between phoenix and 
ATLAS model atmospheres at the F9V+K4V spectral types. 
The results are summarised in Table [3] 
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